Entanglement Entropy for a Coupled Harmonic Oscillator 
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We investigated the entanglement for a simple model of two-particles system in the harmonic 
confinement. The Feynman path integral approach was applied to compute the density matrix of 
the system. The spatial entanglement of ground-state was study using the linear entropy. We found 
the spatial entanglement increases while decreasing of the confining potential uj as this implies that 
the spatial separation between the particles increases. 



PACS numbers: 

I. INTRODUCTION 

Quantum theory is one of the most important branches 
in physics and is well know for being counter-intuitive 
and its bizarre phenomena (nothing like we experience in 
everyday life). Quantum entanglement, which is the di- 
rect consequence of the EPR paradox [J , is a very unique 
phenomenon and has no classical counterparts. Since 
the discovery of quantum entanglement, a wide range 
of applications has been proposed, i.e., quantum cryp- 
tography, quantum teleportation and superdense coding 
[2]. These applications indicate that entanglement could 
be exploited to preform the impossible tasks which can- 
not be done by classical computers. Here is another 
new branch of physics, which is called quantum informa- 
tion and computing, has been born by joining quantum 
physics and theoretical computer science. In this respect 
many solid state systems, especially quantum dots, have 
been proposed as promising hardwares to perform quan- 
tum information processors [3. The main difficulty of 
studying many-body solid state system is the complex- 
ity. The approximation methods must be introduced to 
obtain physical properties, see |4|. 

In this paper, we study the entanglement 0j of the 
ground state of two-body systems using Feynman path 
integral [5]. We will concentrate on the entanglement 
generated by the spatial degrees of freedom and use a cou- 
pled harmonic oscillator as a toy model. The Feynman 
path integral method will be used to compute the den- 
sity matrix of the composite system (bipartite system). 
This provides the way to quantify the entanglement by 
calculating the linear entropy [6]. 

The organisation of this paper is the following. In Sect. 
|ll]we introduce the toy model: a coupled harmonic oscil- 
lator and how to compute the density matrix. Then, the 
linear entropy will be briefly discussed and will be used to 
quantify the spatial entanglement of the system in Sect. 
|III| And the last section will be devoted to summary and 
discussions. 
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II. THE SYSTEM 
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FIG. 1: The system of a coupled harmonic oscillator. 

Consider the coupled harmonic oscillator as show in 
Fig. [1] The particle with mass m is trapped in the har- 
monic potential with frequency w. The particle with 
mass M vibrates around its equilibrium with angular fre- 
quency n. The interaction between these two particles 
is modeled by Hooke's law and k denotes the strength of 
the interaction. The Lagrangian is given by 

^.,y - fi^(r)-|c.V(r) + ^2)^(r)-^f^V(r) 

~\^ir)-y{r)\\ (1) 

where x{t) is the position of the mass m, y(r) is the 
position of the mass M and r is time. The propagator of 
the composite system is then given by 

X,,y = Xixit),y{t),t;xiO),yiO),0) 

= / D(a;)/ 2)(y)e^^ (2) 
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where the action is given by 

S ^ dt£a;,y ■ (3) 

We consider the case of the imaginary time: t — i/3 
where = l/k^T and T is temperature and fee is Boltz- 
mann constant [7]. The propagator Eq. ^ becomes the 
density matrix. 

p^^y = / Vix) / T>{y>-'' , (4) 
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where 



where 
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The reduced density matrix can be obtained by tracing 
out the variable y of mass M [S]. 
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III. LINEAR ENTROPY 

The entanglement generated by the continuous spatial 
degree of freedom is expected to depend on the interplay 
between the strength of the confining potential and in- 
teraction. Then we are interested to investigate how the 
spatial entanglement changes with the system parameter. 

The linear entropy of the reduced density matrix 
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where < < 1 is used to measure the entanglement 
for a pure state and continuous case. For the continuous 
variable case, can be computed by 



Pied = / dxprcd{x,x')prcd{x' ,X) , 
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To obtain Eq. ( 11 1, we take the limit where /3 — )> cxd in 
order to ensure that the system is the ground state. In 
this limit, the system is maximally entangled state when 
§L = 1 leading to the condition that is 
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1 



1 1 

M ^ m 



In the case of the separable state, — 0, we obtain the 
condition 
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Eq. (13 1 and Eq. (14) show that the confining potential 



strength uj (external parameter) depends on the parame- 
ter and K for the case of maximally entangled and the 
case of separable state, respectively. 

Let us, we consider for the case that fl becomes very 
large. We obtain 
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Obviously, the confining parameter uj solely depends on 
the strength of the variable il. On the other hand, if 
becomes infinitesimally small, we obtain 
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and I 



16M 



In this case we see that confining parameter depends 
on the strength of the interaction. From these results 
we can adjust the entanglement between particle mass m 
and particle mass M by tuning the strength of uj with 
vibrating frequency fl and coupling constant k. 



and 
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The system would be completely disentangled where 
Sl ~ 0. On the other hand, the system would be maxi- 
mally entangled state if = 1- 

The linear entropy of our systems is 
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FIG. 2: Linear entropy versus uj 
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FIG. 3: Linear entropy versus fie// 



In Fig.[2]shows the numerical results of how changes 
with respect to w. It is clear the spatial entanglement in- 
creases when the confining potential w decreases. This 
is because the spatial separation between m and M in- 
creases due to low confining potential and hence i^i" 
creases [1]. In Fig. [sj the confining potential uj is fixed 
while the interaction strength k is varied. It is found that 
the spatial entanglement increases with the strength of 
the interaction. 



density matrix of the system could be exactly computed 
by using Feynman path integral approach. We obtain the 
condition on system parameters for maximally entangled 
and completely disentangled state. The system in Sect, 
[njcan be extended by adding more of particles with mass 
M. Then this can be viewed as the particles with AI con- 
stitute the lattice structure while the particle with mass 
m is the charge moving through the lattice. The spatial 
entanglement between the charge (particle with mass m) 
and the lattice (particle with mass M) could possibly be 
studied using the same technique in this paper. Further- 
more, this method could also possible be used to study 
the entanglement in polaron systems. 




IV. CONCLUSION 



We have studied the spatial entanglement through the 
linear entropy for the coupled harmonic oscillator. The 



FIG. 4: The moving particle in the lattice. 
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